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We study the scaling behavior of the step scaling function for SU(3) gauge theory, employing the Iwasaki 
gauge action and the Luescher-Weisz gauge action. In particular, we test the choice of boundary counter terms 
and apply a perturbative procedure for removal of lattice artifacts for the simulation results in the extrapolation 
procedure. We confirm the universality of the step scaling functions at both weak and strong coupling regions. 
We also measure the low energy scale ratio with the Iwasaki action, and confirm its universality. 



1. Introduction 

Recently CP-PACS and JLQCD Collabora- 
tions have started a project for A^f = 3 QCD 
simulations. One of the targets of the project 
is to evaluate the strong coupling constant a-^ 
in the A^f = 3 QCD using the Schrodinger func- 
tional(SF) schemepp. In the project Iwasaki 
gauge action[2] has been employed to avoid the 
strong lattice artifacts of the plaquette gauge ac- 
tion found in A^f = 3 simulations [2. 

In a previous study 0] , as our first step toward 
evaluation of ajj^ for At ^ 3, 0(a) boundary 
improvement coefficients in the SF scheme have 
been determined for the various improved gauge 
actions up to one-loop order in the perturbation 
theory. As the next step, we investigate the lat- 
tice cut off dependence of the step scaling func- 



tion(SSF) non-pcrturbativcly in quenched lattice 
QCD simulations with the improved gauge ac- 
tions for some choices of the improvement coef- 
ficients. We also confirm the universality of SSF 
and the low energy scale ratio, by comparing our 
results with the previous ones obtained by AL- 
PHA CoUaborationiMEl- We refer to ref. |H1 for 
unexplained notations in this report. 

2. Setup 

We follow the same SF setup as in the 
case of the Wilson plaquette action jSj, but em- 
ploy the improved gauge actions including the 
plaquctte(iSo) and rectangle (iSi) loops: 
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The assignment of the weight factor VF,;(C,(?g) 
near the boundary in the time direction is im- 
portant to achieve 0(a) improvement in the SF 
scheme. We employ the weights, 



cocf (.go) = co(l + c. 
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Figure 1. The cut off dependence of E (left-hand side) and Si (right-hand side) at the weak u = 
0.9944(upper) and the strong u = 2.4484(bottom) couphngs for the various improved gauge actions. 
The points ai a/L = in (b) and (d) represent continuum values for the Iwasaki action(open circle), the 
LW action(fulled circle) and also the plaquette actionfcross^ i5l6j . The dotted lines and the straight lines 
represent the combined fit function for the Iwasaki action and the liner fit function for the LW action 
respectively. 



cicf (50^) = ci(3/2 + cf (^^gg + OigD), (3) 

where cq + 8ci — 1 , for the plaquette and rectan- 
gular loops (having 2-link for the spatial direction) 
touching the boundary respectively. The leading 
terms necessary for the tree-level 0(a) improve- 
ment is taken from Ref.[H]. 

The one-loop terms have to satisfy the fol- 
lowing relation to achieve one- loop 0(a) improve- 
ment: cqc^'"^'' -I- 4cicf^*'^'' = Ai/2 where Ai is 
a coefficient of g^a/L term in the SF coupling. 

H(l) P(l) 

We consider two choices: q = 2c^ (condi- 
tion A), and cf^^^ = (condition B). While the 
difference between the conditions A and B is an 



O(a^) contribution at one- loop level, it may be- 
come larger at higher- loop orders. 

Let us introduce quantities which we calculate. 
The SSF describes the evolution of the running 
coupling under a finite scaling factor(say s = 2). 
S denotes the SSF calculated on the lattice. In a 
step to fix the energy scale, one have to set the low 
energy scale ratio imax/'^o a-t the continuum limit 
where Lmax is defined as 5^(imax) — 3-480 and 
To is the Sommer's scale. Later on, we shall show 
the simulation results of their cut off dependence 
and take the continuum limit. 



3. Results 

Simulations for SSF and tq j a are performed on 
CP-PACS. Wc refer to Ref.|Sj for details of simu- 
lation parameters and the continuum extrapola- 
tion procedure. 

In Fig.m we show the cut off dependence of the 
lattice SSF for the improved gauge actions with 
some choices of the boundary counter terms. The 
panels (a) (for weak coupling u — 0.9944) and 

(c) (for strong coupling u = 2.4484) show the 
raw data for E, while a perturbatively improved 
observable SifTD^ is plotted in (b) (for weak) and 

(d) (for strong). The latter is defined as 

Si(2, it,a/L) = S(2,u,a/L)/(l-|-(5i(a/L)u), (4) 

where 5i{a/L) is the one-loop relative 
deviation^. We observe that Si shows a better 
scaling behavior at both coupling regions than 
E. It is particuraly so for the case of the Iwasaki 
action with the tree level 0(a) improvement and 
the condition B. Therefore, we use Si in the 
continuum extrapolation. For the Iwasaki action 
the continuum value is obtained by a joint fit of 
the three sets of data. For the LW action, we 
take the value obtained by a linear fitting of the 
data with the condition A. The continuum values 
are plotted at a/L = in Fig. (b) and (d). We 
observe that the three values obtained with the 
Iwasaki and LW actions in the present work and 
that of ALPHA Collaboration .6^ are consistent 
within la (2.3cr) at the weak (strong) coupling. 

We also show the low energy scale ratio for the 
Iwasaki action with both tree-level and one-loop 
0(a) improvements in Fig. 12 To extrapolate to 
the continuum limit, we use the same fitting form 
as in the case of SSF at strong coupling. We 
do not include the point L/a — 4 for the tree 
level 0(a) improved case. We observe that the 
extrapolated value agrees with the one obtained 
by the plaquette action[2j within errors. 

4. Conclusions 

We have confirmed the universality of the SSF 
and the low energy scale ratio. In the extrapola- 
tion procedure, the perturbative removal of lat- 
tice artifacts well reduces the scaling violation 



1 1 1 




1 1 


1,' 

o 








A' 




- — A' ' 


,■-0 ,. 




1 1 1 


^ - - 


1 1 


m 



0.05 0.1 0.15 0.2 0.25 0.3 

Figure 2. The cut off dependence of the low en- 
ergy scale ratio for the Iwasaki action. The sym- 
bols are the same as the case of the SSF. 

of the SSF for the Iwasaki action with the tree 
level 0(a) improvement and the one- loop 0(a) 
improvement with the condition B. 

As mentioned in the introduction, this work is 
the second step toward unquenched simulations. 
The present study shows that we should use the 
above two choices in the future simulations with 
dynamical quarks. 
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